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ABSTRACT 

Multi-epoch imaging of the Orion equatorial region by the Sloan Digital Sky 
Survey has revealed that significant variability in the blue continuum persists 
into the late-M spectral types, indicating that magnetospheric accretion pro- 
cesses occur below the substellar boundary in the Orion OBI association. We 
investigate the strength of the accretion-related continuum veiling by comparing 
the reddening-invariant colors of the most highly variable stars against those of 
main sequence M dwarfs and evolutionary models. A gradual decrease in the g 
band veiling is seen for the cooler and less massive members, as expected for a 
declining accretion rate with decreasing mass. We also see evidence that the tem- 
perature of the accretion shock decreases in the very low mass regime, reflecting 
a reduction in the energy flux carried by the accretion columns. We find that the 
near-IR excess attributed to circumstellar disk thermal emission drops rapidly 
for spectral types later than M4. This is likely due to the decrease in color con- 
trast between the disk and the cooler stellar photosphere. Since accretion, which 
requires a substantial stellar magnetic field and the presence of a circumstellar 
disk, is inferred for masses down to 0.05 we surmise that brown dwarfs and 
low mass stars share a common mode of formation. 
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1. Introduction 

Brown dwarfs arc defined as objects that have less than the minimum mass required for 
stable Hydrogen fusion. In the absence of rotation and for solar metallicity this mass, referred 
to as the hydrogen-burning limit, is 0.075 Mq or 79 Mj^p (Chabrier & Baraffe 2000). Young 
brown dwarfs are fully convective with photospheric temperatures comparable to M stars. In 
the models of Baraffe et al. (1998) 2 Myr old proto-brown dwarfs of masses 0.020 M© up to 
the hydrogen-burning limit have temperatures between 2500 K and 2900 K, corresponding 
to spectral types of M8.5 to M6 (Luhman 2003). The presence of a circumstellar disk in 
combination with the ability of fully convective objects to generate significant magnetic fields 
should result in young brown dwarfs exhibiting Classical T Tauri behavior (Jayawardhana, 
Mohanty, & Basri 2002). 

The young brown dwarfs in this paper are found as isolated objects within the Orion 
OBI star formation region. Distances and inferred ages of the Orion OBla,b,c, and d 
subassociations are given in Table 1 (from Sherry (2003)). There are two possible scenarios 
for brown dwarf formation: (1) individual formation out of isolated low-mass molecular 
cores, just as stars are presumed to form from higher-massed cores, and (2) ejection of the 
lowest massed objects from a cluster though dynamical interactions. In the latter scenario, 
objects are ejected before accreting enough material to grow into stars, thereby forming 
brown dwarfs (Reipurth & Clarke 2001; Kroupa et al. 2003). 

Consequences of the ejection process include a broader spatial distribution, increased 
velocity dispersion for the lower massed objects, and a truncation of the circumstellar disk 
via tidal forces. The ejection velocity dispersion is of order 5 km/sec corresponding to a 
dispersion in proper motion of 2.6 mas/yr at the 400 pc distance of Orion. Models of the 
ejection process (Armitage & Clarke 1997) predict that the disk will be truncated to an outer 
radius of 10 AU or less, resulting in disk lifetimes of less than 1 Myr. In order to differentiate 
the star-like and ejection formation scenarios we examine observational signatures, such as 
Classical T Tauri behavior, that are due to the presence of a circumstellar disk. 

1.1. Accretion Signatures in T Tauris 

T Tauri variables were originally identified by Joy (1945) as a class of irregular variable 
stars exhibiting marked changes in brightness and color on timescales of hours to days. These 
changes are greatest in the blue and ultraviolet (Herbst et al. 1994), and are expected to be 
related to disk accretion. Magnetospheric accretion models (Konigl 1991) predict that the 
inner circumstellar disk is truncated by the magnetic field of the star. For accretion to occur 
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the truncation must be at or within the Keplerian co-rotation radius, at which the orbital 
period of the disk matches that of the star's rotation. 

Disrupted disk material flows towards the star in accretion columns that follow the 
magnetic field lines connecting the inner edge of the disk with the high latitude regions of 
the star. As evidenced by the broad velocity profiles observed in Ha and other spectral 
features, much of the permitted line emission arises in these columns. Other emission lines, 
including He I and the broad component of the Ca II IR triplet, plus the veiling blue and 
ultraviolet continua arise in the high temperature shock fronts at the base of the accretion 
cohimns near the stellar surface (Gullbring et al. 2000). The lifetime of the classical T Tauri 
phase is thought to be 1 to 10 Myr (Kenyon & Hartmann 1995). 

Protostellar magnetospheric accretion is one of a dynamic process characterized by 
intricate magnetic field topologies and unstable field-disk interactions. Detailed MHD com- 
putations by Romanova (2002), Kiiker, Henning, & Riidiger (2003), and von Rekowski & 
Brandenburg (2003) indicate that even for the ideahzed case of a dipole field aligned with the 
stellar rotation axis the accretion process is fundamentally unstable. An intensive photomet- 
ric and spectroscopic campaign targeting the nearly edge-on Classical T Tauri AA Tau by 
Bouvier (2003) shows clear evidence for large scale instabihties developing in T Tauri magne- 
tospheres as the magnetic field lines are twisted by differential rotation between the star and 
the inner disk. Observational support of magnetospheric accretion in AA Tau includes time 
delays between the Ha, H/?, He II line emission and the accretion shock generated emission 
consistent with free-fall from ~8 stellar radii and the presence of two rotationally modulated 
hot spots. 

1.2. Accretion in Very Low Mass Stars and Brown Dwarfs 

There has been considerable effort over recent years to observe and characterize the 
formation processes in young very low mass stars and brown stars. This has included analyses 
of spectroscopic signatures of accretion, veiling measurements, variability studies, and IR 
searches for circum(sub)stellar disk emission. 

1.2.1. Emission Lines 

Ha emission lines with widths in excess of 200 km/s have been detected in young brown 
dwarfs having masses as low as ~ 20 Mj^p in many star formation regions including IC 348, 
Taurus, Chamaeleon I, Upper Scorpius, and p Oph (Jayawardhana, Mohanty, & Basri 2002; 
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MuzeroUe et al. 2003; Natta et al. 2004; Mohanty, Jayawardhana, & Basri 2005; MuzeroUe 
et al. 2005). The mass accretion rate appears to scale roughly as which is steeper than 
the linear relation inferred by White & Basri (2003) based on studies in Taurus-Auriga. In 
addition, the uncertainties in the determination of stellar parameters (mass, radius) from 
evolutionary models of very low mass objects (Baraffc et al. 2002) during first few Myr 
contribute another source of uncertainty in these relations. However, Padoan et al. (2005) 
propose that the accretion rate is controlled by accretion from the surrounding gas in the star 
formation region rather than from the disk evolution. The resulting Bondi-Hoyle accretion 
model yields trends consistent with the MuzeroUe et al. (2003) relation although the vahdity 
of this hypothesis needs to be checked through detailed observations and simulations. 

Natta et al. (2004) find that the mass accretion rate scales exponentially with the width 
of the Ha line measured at the 10% intensity level by comparing previous results in the 
substellar domain (MuzeroUe et al. (2003), White & Basri (2003), and Barrado y Navascues 
et al. (2004)) with observations in Chamaeleon I and p Oph and study of Classical T Tauri 
stars(Gullbring et al. 1998). The lowest detected accretion rates are on the order of 10~^^ 
Mo/yr. 

1.2.2. Optical Veiling 

Optical veiling measurements of young brown dwarfs have been primarily obtained red- 
ward of 6000 A. Barrado y Navascues, Mohanty, & Jayawardhana (2004) find r62oo = 1-0 and 
^6750 — 0.25 from the infilling of several TiO bands in the very low mass object LS-RCrA 
1, where is the ratio of the excess emission to the photospheric emission measured at a 
wavelength A (A). Due to the presence of [S 11], [O I], [O II], and [N II] forbidden lines in 
the spectra they propose that this substellar object is driving an outflow analogous to that 
which forms the Herbig-Haro shocks seen near intermediate and low mass T Tauri stars. 

Accretion shock models for very low mass stars (MuzeroUe ct al. 2000) suggest that 
accretion rates as low as 10~^ MQ/yr would be detectable in the Johnson V band against 
an M6 photosphere. We expect a greater sensitivity to accretion rates when working in the 
Sloan Digital Sky Survey (SDSS) g band due to increased brightness contrast between the 
shock and the photosphere. 
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1.2.3. Variability 

Optical photometric variability of low mass pre-main sequence stars is the result of four 
physical effects. These are rotational modulation of cool spots, rotational modulation of 
hot spots formed at the base of magnetospheric accretion columns, instabilities in the mass 
accretion rate, and flaring. 

The study by Herbst et al. (1994) finds three distinct patterns of variability in Weak- 
lined T Tauris (WTTS) and Classical T Tauris (CTTS). The Type I variables are charac- 
terized by low amphtude multi-band fluctuations on timescales of ~ 0.5 to 18 days. This 
behavior is primarily seen in WTTS but can been detected in some WTTS. The second type 
of variation only occurs in the CTTS and is marked by ireegular high amplitude brightness 
changes, especially in the near ultraviolet and the blue, on timescales as short as a few hours. 
The Type III variations, like the Type II, are irregular but occur on much longer timescales 
of days to weeks. This last category are only seen in pre-main sequence stars having early 
spectral types (AO to Kl). 

Type I variations are believed due to rotational modulation of cool spots on the stellar 
surface. These spots are a consequence of the magnetic activity and can persist for over 100 
rotational periods yielding periodic variations in VRI. Flaring is often seen in the B and 
especially U bands. WTTS only exhibit Type I variations. When magnetospheric accretion 
is present, i.e. in the CTTS stage, then Type II variations are seen. These are due to a 
combination of rotational modulation of the accretion hot spots together with short period 
fluctuations in the mass accretion rate. In most cases the variations in the accretion mask 
the rotation signature. Herbst et al. (1994) also identify a third class (Type III) variables 
such as RY Tau and SU Aur that exhibit large non-periodic variation in the V band but do 
not show significant veiling. This behavior is thought to be caused by variable obscuration 
but is also limited to intermediate and high- mass stars (spectral types AO to Kl); we do not 
expect to find Type III variables in our low-mass sample. 

In this sample of intermediate and low mass T Tauris Herbst et al. (1994) find that the 
maximum amplitudes seen in Type I variables are 0.8 magnitudes in V and 0.5 magnitudes 
in I for V410 Tau (SpT = K3). For the Type I variables of spectral type MO and later the 
V band amplitudes are typically less than 0.3 magnitudes. 

Variability studies of young very low mass stars and brown dwarfs have primarily used 
the longer wavelength filters, most notably the Johnson-Cousins /. Herbst et al. (2002) 
utilized an intermediate band filter centered at 815.9 nm to study variability in the Orion 
Nebula Cluster and found that 46% of the stars with a < 0.1 were periodic, where a is the 
standard deviation measured over all observations for a star. Only 24% having cr > 0.1 were 
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periodic, the majority exhibited the irregular variation characteristic of Type II variables. 

In their study of IC 348 Cohen, Herbst, & WiUiams (2004) found that the full range 
of variation for periodic stars of spectral type M0-M4 or later is < 0.6 magnitudes in the / 
band. No CTTS exhibited periodic variability over 5 seasons of monitoring. Cohen, Herbst, 
& WiUiams (2004) also characterized the variabihty in their sample by a finding that in their 
/ < 14.3 subsample the WTTS had a < 0.1 and the CTTS were found to have a > 0.05. 
Only the Type II variables (CTTS) had a > 0.1. 

The survey of the Cha I star formation region in R, i, and J by Joergens et al. (2003) 
found that the largest amplitudes in each band were 0.18, 0.14, and 0.13, respectively. The 
earliest spectral type in the sample was M5, corresponding to a maximum mass of 0.12 Mq. 
These objects showed periodic variations, presumably due to spots, whose peak-to-peak 
amplitudes drop in later spectral types. 

I band monitoring of the a Ori and e Ori clusters were performed by Scholz & Eisloffel 
(2004a) and Scholz & Eisloffel (2005). They found that periodic variables had amplitudes 
less than 0.2 magnitudes. The larger amplitude variables were generally non-periodic. S Ori 
45, a ~ 20 Jupiter mass member of the a Ori cluster (SpT = MS. 5) was studied by Zapatero 
Osorio et al. (2003) in the IJ bands who measured a periodic peak to peak amplitude ~ 0.2 
magnitudes. 

Caballcro ct al. (2004) find that nine out of 32 young brown dwarf candidates (spectral 
types M5.5 to L2) exhibit / band variability. The amplitudes of mid-term (day-to-day) varia- 
tions were as high as 0.36 magnitudes while shorter term variations had standard deviations 
~ 0.05 magnitudes. Correlations were found between high amplitude variability (> 0.12 
magnitudes) and detection of either a near-IR excess or strong Ho; emission indicative of a 
possible accretion disk. 

1.24. Near-IR Excess 

The thermal emission from the inner circumstellar disk can be detected in the near-IR 
H and K bands (Meyer, Calvet, & Hillenbrand 1997). For the less massive and cooler young 
brown dwarfs the reduced contrast between the disk emission and the substellar photosphere 
motivates disk surveys in L band. An L band survey by Liu, Najita, & Tokunaga (2003) 
finds that for low mass stars and substellar objects in IC 348 and Tauriis the disk fraction 
is not dependent on mass or age and that the presence of disk emission is associated with 
accretion indicators such as Ha emission. 
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Comparison of disk models with mid-IR photometry shows that the disks surrounding 
young brown dwarfs can be flared like those around Classical T Tauri stars (Mohanty et al. 
2004). The model fits also suggest the presence of inner holes that are a few substellar radii 
in size. 



1.3. The Orion Equatorial Region 

The Orion OBI a,b,c and d sub-associations have been the targets of several surveys. 
These include the optical photometric and spectroscopic study of the a Ori cluster in Orion 
OBlb by Walter, Wolk, & Sherry (1998), the Ha Kiso object prism survey (Wiramihardja 
et al. 1989), and the ROSAT All-Sky Survey (Alcala et al. 1996). With the exception of the 
Orion Nebula Cluster (OBld) these arc found along the celestial equator and are considered 
in this study. The estimated ages span 1.7 Myr to 11.4 Myr with distances between 330 and 
460 parsecs corresponding to distance moduli of 7.6 to 8.3 magnitudes. 

Recent imaging surveys of the Orion equatorial region include Wolk (1996), Sherry 
(2003), and Briceho (2003). Wolk (1996) performed UBVRcIc photometric monitoring of X- 
ray selected stars in the Orion OBla and OBlb associations with an estimated completeness 
hmit ol V — 18.5. Sherry (2003) conducted a deeper BVRcIc survey covering 5 square 
degrees in Orion OBlb around a Ori, e Ori, and S Ori using the CTIO 0.9 and 1.5 meter 
telescopes. Their adopted faint hmit was V — 20.5. The B band photometry was not 
analyzed. 

The CIDA-QUEST variabihty survey (Briceiio 2003) is acquiring multi-epoch drift- 
scanned imaging in the Orion equatorial region using the 1.0/1.5 meter Schmidt telescope at 
The National Astronomical Observatory of Venezuela. Exposures are acquired in the VRIHa 
bands at a plate scale of 1.02" per pixel. The lOo" limit is approximately V = 19.7 although 
in selection of candidate pre-main sequence stars they specify V — 16 — 18.5 (Briceno 2004), 
resulting in a minimum mass limit of ~ 0.12 Mq. 

Studies of neighboring regions include Barrado y Navascues et al. (2004) who performed 
a deep Rcic survey of the A Orionis cluster (d = 400 pc; 5 Myr) with the CFHT 12K camera. 
Their completeness limits were 22.75 in both bands with the lowest mass confirmed cluster 
member having 0.02 Mq with a spectral type of M8.5. Studies of NGC 2264 including S 
Mon and the Cone Nebula region (RebuU et al. 2002) and of the Orion Nebula Flanking 
fields (RebuU et al. 2000) explicitly used ultraviolet and near-IR excess as disk indicators. 
Their hmiting magnitudes are [/ ~ 20, F ~ 20, and Ic ~ 18. 
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1.4. Program 

The goal of this paper is to study the trends in veihng emission, variabihty, and near- 
infrared excess in the low mass stars and brown dwarfs found in the Orion OBla,b, and c 
subassociations. If brown dwarfs form in the same manner as stars we expect that these 
characteristics will vary smoothly across the stellar /substellar mass boundary. 

The remainder of this paper is organized as follows: the SDSS and Two Micron All-Sky 
Survey (2MASS) observations used in this paper are presented in §2. In §3 we discuss the 
expected colors of young mass stars and brown dwarfs using the colors of spectroscopically 
verified field M dwarfs (West et al. 2004), the evolutionary models of Baraffe et al. (1998) 
(hereafter BCAH98), and the bolometric corrections of Girardi et al. (2004). The selection 
criteria for young low mass objects including an estimate on the contamination fraction are 
presented in §4. We discuss the empirical trends in variability and color in §5 and in §6 we 
provide a summary of results and upcoming work. 

2. Observations 
2.1. Photometry 

The Sloan Digital Sky Survey (SDSS) obtains deep photometry with magnitude limits 
(defined by 95% detection repeatability for point sources) of -u = 22.0, g = 22.2, r = 22.2, 
i — 21.2 and z = 20.5. These five passbands, ugriz, have effective wavelengths of 3540, 
4760, 6290, 7690, and 9250 A, respectively. A technical summary of the SDSS is given by 
York et al. (2000). The SDSS imaging camera is described by Gunn et al. (1998). Ivezic et 
al. (2004) discuss the data management and photometric quality assessment system. 

The Early Data Release and the Data Release One are described by Stoughton et al. 
(2002) and Abazajian et al. (2003). The former includes an extensive discussion of the data 
outputs and software. Pier et al. (2003) describe the astrometric calibration of the survey 
and the network of primary photometric standard stars is described by Smith et al. (2002). 
The photometric system itself is defined by Fukugita et al. (1996), and the system which 
monitors the site photometricity by Hogg et al. (2001). Abazajian ct al. (2003) discuss 
the differences between the native SDSS 2.5m ugriz system and the u'g'r'i'z' standard star 
system defined on the USNO 1.0 m (Smith et al. 2002). 

The SDSS low Galactic latitude data which includes the Orion equatorial imaging used 
in this work are described by Finkbeiner et al. (2004). 2MASS (Skrutskie et al. 1997) 
obtained nearly complete coverage of the sky in JHKg. The multi-epoch data (Table 2) we 
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use in this study were obtained under photometric conditions and cover 80° < q;2ooo < 90° 
and -1.25° < ^onno < 1-25°. 



2.2. Reddening-Invariant Indices 

In this work we employ reddening-invariant indices of the form Qxyz = — y) — {y — 
z) X E{x — y)/E{y — z). Qxyz is dependent upon the assumed ratio of general to selective 
extinction {Ry = Av/E{B - V); CardeUi, Clayton, & Mathis (1989)). Here xyz refer to 
the specific passbands, e.g. ugrizJHKg and E{x — y) is the color excess due to reddening 
in the x — y color. This definition of reddening-invariant colors follows the original Johnson 
& Morgan (1953) Q that would be written as Qubv in our notation. In the {x — y,y — z) 
color-color diagram the Qxyz axis is perpendicular to the reddening vector. 

This approach has previously been used to study prc-main sequence populations by de 
Geus, de Zeeuw, & Lub (1989) and Brown, dc Gcus, & dc Zccuw (1994) for the Scorpius- 
Centaurus OB and Orion OB lb associations, respectively. In both cases reddening-invariant 
colors were formed by [x — y] = {x — y) — {V — B) x E{x — y) / E{V — B) . using the VBLUW 
Walravcn system. Stellar parameters were estimated using reddening-invariant two-color 
diagrams to compare observed stellar colors against those from Kurucz model grid (Kurucz 
1979). 

While the contrast between the accretion shock and the cool photosphere is greatest 
in the SDSS u band, we choose not to use the Qugr color due to sensitivity issues and the 
presence of a "red leak" in the SDSS imager u filter. M dwarfs are 2.5 magnitudes fainter 
in u than in g, thus use of the u band in Orion would not reach young brown dwarfs. 

We use the extinction tables derived by D. Finkbeiner^ to define the coefficients used in 
defining reddening-invariant colors. These tables contain the Ax/E{B — V) values for the 
SDSS ugriz filters for specific values of Ry and source spectra. The values we present here 
are obtained using an F dwarf source spectrum and Ry — 3.1 and 5.5. The Ry — 3.1 case 
is the standard extinction law found in the diffuse ISM. The Ry — 5.5 law is representative 
of that found in molecular clouds due to larger dust grains and is shown here for example. 
Throughout this paper we will adopt Ry = 3.1 as it is characteristic of all but perhaps the 
most heavily extincted regions of the Orion OBlb association (Warren & Hesser 1978). 



private communication; see http://www.astro.princeton.edu/~dfink/sdssfilters/ 
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\{9-r)- 1.339(r - 1) Ry = 5.5 

and 

^ r (r - i) - 0.987(i -z) Rv^ 3.1, 



The near-IR reddening- invariant color is independent of the value of Ry, from Schlegel, 
Finkbeiner, & Davis (1998) we obtain 

QjHK = {J-H)- l.5m{H - K). (3) 



3. Expected Colors of Young Brown Dwarfs 

In this section we compare the expected colors of young brown dwarfs against those 
of field M dwarfs of equivalent temperatures. Examination of the temperature - surface 
gravity (Figure 1) relations based on the BCAH98 models shows for objects at the Hydrogen 
Burning Limit younger than 10 Myr the effective temperature ranges between 2900 K and 
3000 K and the surface gravity varies from log(g) = 3.5 to 4.2. 

When the locations in Hertzsprung-Russell diagrams of spectroscopically classified and 
presumably co-eval stars in specific star formation regions are compared against evolutionary 
models the inferred ages mismatch using the spectral type to Tg// scale defined for disk 
dwarfs. This has led to the adoption of a semi-empirical temperature scale intermediate 
between that of the dwarfs and the giants (Luhman 2003). This practice is motivated by the 
surface gravities expected in pre-main sequence stars and results in effective temperatures 
several 100 K warmer at a given spectral type. Thus a T^ff of 2900 K corresponds to a 
spectral type of M4.5 for a dwarf (Reid & Hawley 2000) but of M6.5 for a pre-main sequence 
star. 

The photospheric colors of M dwarfs and M giants differ at the later spectral types 
with the M giants becoming bhier due to surface gravity effects. This color shift with lower 
surface gravity is due to the change in the optical depth for the formation of molecules, 
most importantly TiO, that dominate the stellar spectrum. Luhman (1999) notes that for 
objects of later spectral types, e.g. M8, the trend is reversed with the lower surface gravity 
stars becoming redder longward of 8500 A. Due to these deviations from disk dwarf colors 
UV-excess and blue-excess techniques based on the ZAMS become ambiguous for pre-main 
sequence stars later than M4 (RebuU et al. 2000). 



3.1. Disk Dwarf Colors 



Figure 2, based on the sample of West et al. (2004), shows that the Qriz color index 
increases towards cooler spectral types among M dwarfs, reaching a maximum value ~ 1.2 
at M8. At the ages in which we expect to see accretion activity this turnover occurs for 
objects near 0.02 Mq which are too faint to be included in this survey. 

To study the change in Qgn relative to the disk dwarf locus we obtain linear fits be- 
tween Qgri and Qriz using photometry of M dwarfs in the SDSS Third Data Release (DR3; 
Abazajian et al. (2005)) obtaining 

^ r-3.371Q„, + 1.119 i?v = 3.1, 
^gri I -2.24:lQriz + 1.163 Ry = 5.5 ^ ^ 

using a minimum threshold to Q^z of 0.5 to avoid the slope change in the color-color diagram 
at the earliest M dwarfs. The fit residuals have RMS values of 0.148 and 0.117 for Ry — 3.1 
and 5.5, respectively. 

The Rv reddening-invariant two-color diagram is shown in Figure 3 with the theoretical 
isochrones of Girardi et al. (2004) overplotted. These isochrones are based on the "AMES" 
(4000 K > Te// > 2800 K) and "AMES-dusty" (2800 K > Te// > 500 K) model atmospheres 
of Allard, Hauschildt, & Schwenke (2000) for cool stars and brown dwarfs. While the latter 
cover log(g) between 3.5 and 6.0, the ugriz bolometric corrections were only computed for 
the warmer "AMES" models having log(g) of 5.5 and 6.0. 

The fit residuals (AQ^^j) are plotted against Qriz for the model isochrones and the 
spectroscopic sample of West et al. (2004). For the purposes of illustration this M dwarf 
sample is subdivided at M4.5 with the two resulting histograms in color-color space scaled 
to match peak counts. As is seen in Figure 4 the observed reddening-invariant colors are 
well-described by the log(g) = 5.0 and 5.5 models appropriate for the disk dwarfs. The 
blueward trends in both Qriz and AQgri increasing with decreasing surface gravity as shown 
by the log(g) = 4.5, 4.0, and 3.5 models. The dotted lines connect models of Te// — 2500 
K, 2700 K, 2900 K, and 3100 K. We refer to AQgri as the observed Qgri offset since it is 
measured relative to the {Qriz, Qgri) linear fit. As described below, understanding the effect 
of low surface gravity on AQgri is crucial for estimating the g band veiling from reddening- 
invariant colors. 
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3.2. Veiling Models 

The identification of accretion shocks by photometry involves detection of the additional 
continuum fiux generated by the reprocessing of the high temperature shock back along the 
accretion column. The veiling produced by this excess continuum is defined as the ratio of 
the shock flux to that of the stellar photosphere within a specifled wavelength range or filter 
passband. 

The observed increase in brightness in an arbitrary filter x is 

Am, = -2.51og(l + r,) (5) 

where Tx is the veiling in x and nix is the observed magnitude. The observed veilings in low 
mass T Tauris are significantly higher in the near- ultraviolet and blue regions of the spectrum 
than in the red. For all but the mostly heavily veiled sources the contribution of the accretion 
continuum to the photosphere is minimal at wavelengths above 7000A (MuzeroUe et al. 2003) 
and therefore should be minimal in the SDSS i and z bands. This is a contrast effect between 
the ~ lO'^ K shock emission and the < 3000 K star. 

The change in a reddening-invariant color due to veiling is then 

^Qxyz = -2.5{log[(l+r,)/(l + r,)] + 

log[(l + r,)/(l + r,)]E{x - y)/E{y - z)}. (6) 

Our veiling models are characterized by the veilings in the SDSS g and r band, and r^. 
We explicitly assume that rj and are zero (c/. supra) and obtain 

AQgri = -2.5{log[(l + rg)/{l + r,)] - log[l + rr]E{g - r)/E{r - i)} (7) 

and 

AQH^ = -2.5{log[l + r,]}. (8) 

In the general case where > we can not deduce the photospheric spectral type from 
Qriz- However, since veihng makes stars look bluer (earlier, more massive), we can use the 
observed Q^z to deduce the earliest possible spectral type. For veiling models consisting of a 
non-zero rg and we can examine the variations in the Q^z and Qgri color. Figure 5 shows 
the color shifts for Vg varying continuously from -0.5 to 1 and r^/r^ = 0.00, 0.25, 0.50, 0.75, 
and 1.00. The negative veihng values are used to model the cool spots seen in Weak-lined T 
Tauris. 

When an unknown veiling is present in the r band, the key diagnostic feature in the 
{Qriz, Qgri) color-color diagram is the vertical displacement of the star from the stellar locus 
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along the Qgri reddening-invariant index. This is shown in Figure 6 as a function of Vg 
assuming the hnear fits to the locus given above. In the case of weak to moderate veihng, 
< 1, we see that differing the to Vg ratio only results in a minimal spread of derived 
Tg for a given AQgri. We can therefore use the observed AQg^i as a proxy for continuum 
veiling even when the spectral type and Qriz color are not known a priori assuming that the 
intrinsic Qgri offset due to low surface gravity can be inferred. 

4. Candidate Selection 

Possible low mass and very low mass accreting stars were selected on the basis of color, 
location in color-magnitude diagrams, and variability. An estimate of contamination was 
obtained by applying these criteria to a neighboring field. 

4.1. Color and Magnitude Criteria 

Wc adopted selection criteria (r — i) > 0.6 and Qriz > 0.35 to target M dwarf colored 
stars. This corresponds to T^ff < 3700 K or spectral types later than MO. 5 in the field and 
Ml on the Luhman (2003) scale. The maximum stellar mass inferred from the 2 Myr Baraffe 
et al. (1998) isochrone is 0.8 Mq. 

To reduce the scatter in the reddening invariant colors we required photometric errors 
less than 0.05 magnitudes in grizJH. This was relaxed for the Ks band to extend the 
detection limit. The resulting magnitude and error limits are summarized in Table 3. The 
limiting J band magnitude of 15.5 corresponds to a lower mass limit of 0.05 M© (~ M7) 
based on the BCAH98 2 Myr isochrone. Propagation of these errors into the expressions 
for the reddening-invariant colors yields maximum errors for Qgri of 0.088 and 0.075 and for 
Qriz of 0.045 and 0.046, each in the cases of Ry — 3.1 and 5.5, respectively. The maximum 
error in Qjhk is 0.159 which is dominated by the Kg error of 0.1 at the faint limit. 

In order to remove background stars we required that candidates were brighter than 
the main sequence placed at the distance of Orion in both the (r — i,r) amd {i — z,i) 
color-magnitude diagrams. The color- absolute magnitude relations in the SDSS passbands 
for late K and M dwarfs are Mr = (6.110 ± 0.418) + (3.800 ± 0.320) (r - i) and = 
(5.936 ± 0.496) + (6.308 ± 0.803) (i - z) (Golimowski et al. 2005 in prep.). The errors in 
these relations are dominated not by photometric uncertainly (which is typically 1-2%), but 
by the intrinsic cosmic scatter of the sample. 
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4.2. Variability 

Due to the non-periodic variations of low mass CTTS and the irregular sampling by the 
SDSS we expect the resulting lightcurves to be fairly random. A truly random sequence has 
a standard deviation (cr) equal to 1.0/\/T2 of the peak-to-peak amplitude. Given a finite 
number of samples and a probability that the true value of a exceeds a specific threshold 
we need to find the corresponding minimum threshold for the observed a that indicates 
significant variabihty. We generated a series of Monte Carlo simulations to determine the 
maximum cr in a random sequence given an intrinsic a of the parent population and a limited 
number of observations. The scaling factors presented in Table 4 refiect the maximum a seen 
in 99% of the simulated data sets using 10,000 realizations per point. In the limit of large 
number of observations the scaling factors tend towards unity. 

We identify variable stars on the basis of the g band a exceeding both a threshold of 0.05 
magnitudes (scaled by the correction factor defined above), and three times the photometric 
error in g. Of these we also select a subset that meet the same conditions applied to the z 
band. This z band variability criterion matches the typical boundary between Type I and 
Type II / band variations in very low mass stars and young brown dwarfs. 

4.3. Results 

On the basis of these criteria we identify 507 stars that are significantly variable in the 
g band of which 215 also meet the criteria for z band variability. In Figures 7 and 8 we see 
that the majority of the variable stars form a locus parallel to and roughly 2 magnitudes 
above the main sequence placed at the m — M = 8.2 distance of the Orion star formation 
region. The spatial distribution of these low mass candidates is concentrated in the Orion 
OBla and OBlb associations and the NGC 2068/2071 protocluster (Figure 9). The Orion 
OBlc subassociation members comprise a dispersed population in this region. 

4.4. Compctrison Field 

In order to assess the possible contamination of our Orion sample by variable field stars 
we apply the criteria defined above to a neighboring equatorial field bounded by 40° < a < 
70° and —1.25° < 5 < 1.25° which encompasses three times the area of our survey. This is a 
complex region containing at least three populations of young stars which we now describe. 

On the basis of ROSAT All-Sky Survey (Triimpler 1983) detections Neuhaiiser (1997) 
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subsequently identified a population of ~ 30 Myr low mass WTTS south of the Taurus- Auriga 
cloud. While we expect to see Class I variability in these stars signs of active accretion or 
the presence of inner circumstellar disks should not be evident due to the typical 10 Myr 
disk lifetime. We also expect to detect the background 8 Myr old Gould Belt stars (Guillout 
et al. 1998) which may include CTTS. As the Gould Belt midplane crosses our equatorial 
survey area near a = 70 we suspect some Gould Belt stars might be found in our Orion 
sample. 

The third star formation complex is that associated with the nearby (80 ± 20 pc) 
translucent cloud MBM 18 (Larson & Rccd 2003). This high latitude cloud contains compact 
regions of CO(l-O) emission with line intensities above 3 K. If MBM 18 is an active star 
formation region then due to its proximity we expect all but the lowest massed young stars 
to be too bright for the SDSS, 

Application of our Orion CTTS criteria results in 133 g band variables of which 19 
are also variable in the z band. After scaling by the ratio of the survey areas these imply 
possible contamination fractions of 8.7% and 2.9% for the Orion g band variables and its z 
band variable subset, respectively. 



5. Empirical Results For The Orion Population 

5.1. Veiling Signatures 

In Figure 10 we show the reddening-invariant Qriz color shift between the faint and 
the bright states defined by the g band brightness against the 2MASS J magnitude. The 
observed trend is for Qriz to decrease when the star is brighter in g which is consistent 
with the presence of veiling in the r band. In this figure the horizontal dotted line marks 
the substellar boundary based on the BCAH98 2 Myr isochrone placed at m-M = 8.22. 
We convert the BCAH98 JHK magnitudes from the CIT to the 2MASS system following 
Carpenter (2001). The derived effective temperature from Qriz changes by approximately 100 
K for every 0.1 change in the color (see Figure 2). Due to this uncertainty in determination 
of the photospheric Q^z color we will present all color and variability trends against the 
2MASS J magnitude which should be relatively unaffected by both shock emission and disk 
thermal emission. 

In Figure 11 we see that the bright state AQg^i during both faint and bright states 
slightly decreases in the less massive (fainter) objects. However, when the observed AQg^j 
are compared against those predicted from the cooler (~ 2700 K) low surface gravity models 
(dot-dash hue) it is clear that the bulk of veihng signature is probably due to the shifts 
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in photospheric colors. Thus, at the substellar boundary we only see evidence for weak g 
band veiling which consistent with the results of MuzeroUe et al. (2003) who find a deficit of 
continuum veihng at 5500 A for accreting stars of spectral types M5 to M7. 



5.2. Variability 

If the veihng in a specific band x ranges from r^""' to r^^^^ then the observed magnitude 
change is 



X 



Am^ ^ -2.5109 [j^). (9) 



This scales montonically, albeit non-linear ly, with r^*^'*, if 

^^high - ^ ^ ^high y ) 

which in the case where r^°^ scales with r^*s^^ i_g_ j-iow _ (y^^high Ycqmrcs o; < 1, which is 
implied by r^°'^ < r^*^'*. In general we then expect the amplitude and RMS of the variability 
to be related to the veiling and hence with the shock flux in that band. 

In Figure 12 we see that the maximum RMS of the g band variations peaks at J ~ 13 
and that no high-amplitude (c^ > 0.5) variables are seen below J = 14. The triangles 
mark the median ag values computed in for one-magnitude bins spanning 11 < J < 16. We 
also see that the maximum ratio of the g band to z band RMS {ag/az) peaks at the same 
magnitude and then decreases in the fainter stars (Figure 13). Only five faint (J > 14) 
candidate accreting stars have Ogjuz > 3. This trend is perhaps indicative of cooling of the 
accretion shock temperature with diminished mass. 

For the intermediate mass Classical T Tauris studied by GuUbring et al. (1998), the 
Balmer jump diminishes as the veiling and shock bolometric luminosity increase. Conversely, 
for the weakly accreting systems expected in very low mass stars and young brown dwarfs 
the Balmer jump should be relatively large suggesting that the bulk of the shock emission 
will be in the near-ultraviolet rather than the optical bands. The shock models of Muzcrolle 
et al. (2003) suggest that shocks with lower energy flux carried by the accreting column emit 
more at redder wavelengths due to quasi-blackbody emission from the heated photosphere. 
The trend we see is consistent with a decrease in the shock energy flux with decreasing mass. 
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5.3. Near-Infrcired Excess 

Near-Infrared (JHK) photometry has often been used to detect the presence of excess 
emission primarily in the K band due to thermal emission from the inner circumstellar disk. 
Meyer, Calvet, & Hillenbrand (1997) find a Classical T Tauri locus in {J — H, H — K) based 
on the dereddened colors of K7/M0 intermediate massed stars in Taurus. In Figure 14 we 
see that the Kg excess indicating either a disk signatures or significant extinction drops for 
objects having J < 14. For the 2 Myr BCAH98 isochrone this corresponds to 0.11 M© and 
Te// = 3000 K (~ M5.5). This lack of signature in JHK is expected for substellar objects 
since they are less luminous and so the inner circumstellar disk will be cooler. As discussed 
by Liu, Najita, & Tokunaga (2003) this contrast problem between the disk and the stellar 
photosphere motivates searching for young brown dwarf disks in the L and longer wavelength 
bands. 



6. Summary 

Analysis of 2MASS and multi-epoch SDSS photometry of low mass pre-main sequence 
objects stars within the Orion OBI assocation shows that: 

1. g band variations having a standard deviation greater than 0.05 magnitudes persist 
for objects having reddening- invariant colors redder than, and J magnitudes fainter than, 
that expected for the stellar/substellar boundary. 

2. While identification of g band veihng for spectral types later than M4 is hindered by 
our uncertain knowledge of surface gravity effects on photospheric colors, we see evidence 
for accretion signatures even in these very low mass objects. The general trend is for the 
veiling signature to diminish with decreasing temperature (and mass). 

3. No disk signatures based on J — K color excess are seen for spectral types later than 
M4. We infer based on variability and weak veiling that the accretion process and hence 
the disks are present, and that relative contribution to the K band flux is weak. This is 
consistent with earlier work on young brown dwarf disks. 

We have demonstrated that multi-epoch and multi-band imaging surveys enable identi- 
fication and classification of pre-main sequence objects based on variability and reddening- 
invariant colors. While use of the SDSS g band was mandated by the sensitivity issues and 
the red-leak in the imager's filter system, inclusion of a broad-band filter that is shortward 
of the Balmer jump such the SDSS u in the next generation of surveys, i.e. the Large Syn- 
optic Survey Telescope or LSST, will increase the sensitivity to the veiling emission from the 
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accretion shock for the later and less massive objects. 

We are supplementing this imaging survey with an SDSS spectroscopic program target- 
ing low mass stars in Orion and the Taurus- Auriga star formation region. Results from this 
study in addition to analysis of the veiling trends presented here and their constraints on 
shock luminosities and mass accretion rates will be discussed in subsequent work. 
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Teff (K) 

Fig. 1. — Temperature and surface gravity for young low mass stars. The predicted 
effective temperatures (K) and surface gravities (log(g)) are shown here for the BCAH98 
models having 0.02Mq < M < I.OMq and ages between 1 Myr and 10 Myr. Isochrones 
are drawn for log(years) = 6.0 to 7.0, spaced by 0.2 dex. The thick black line traces the 
Hydrogen Burning Limit (O.OTSM©). 
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Fig. 2. — Model Qriz color and T^ff relations. This figure shows the predicted Qriz 
colors as a function of effective temperatures using the ugriz bolometric correction tables of 
Girardi et al. (2004) which are based on the AUard, Hauschildt, & Schwenke (2000) model 
atmospheres. The lower surface gravity models (log(g) = 3.5, 4.0 and 4.5) were not presented 
by Girardi et al. (2004) for Tg// > 2700. The vertical line marks the approximate location 
of the Hydrogen Burning Limit for the first 10 Myr. 
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Fig. 3. — Disk dwarf reddening- invariant color-color relations. This figure shows 
{Qriz: Qgri) color-color diagrams for 1,116,184 M star candidates in the SDSS Third Data 
Release (Abazajian et al. 2005). The dotted diagonal hne traces hnear fits based on the 
Qriz > 0.5 subsample of 694,062 stars. The Girardi ct al. (2004) grid is overlaid with dashed 
hnes connecting models having Te// = 2500 K, 2700 K, 2900 K, and 3100 K. 
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Fig. 4. — Fit residuals for SDSS M dwarfs. The fit residuals (AQgri) are plotted against 
Qriz for the Girardi et al. (2004) model isochrones and the SDSS spectroscopic sample of 
West et al. (2004). The meaning of the lines follows Figure 3. For the purposes of illustration 
this M dwarf sample is subdivided at M4.5 with the resulting 2-D histograms scaled to match 
peak counts. 
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Fig. 5. — Color changes due to veiling. The changes in Qriz and Qgri are shown for the 
ratio of r,. to Vg varying from 0.0 to 1.0 in steps of 0.25 assuming a Ry = 3.1 reddening law. 
The g band veiling varies from -0.5 to 1.0 along each curve. The dotted line shows the linear 
fit to the disk dwarf locus assumed in computation of AQ^ri- 
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Fig. 6. — Relation between observed AQgri and veiling. The observed vertical separa- 
tion (AQgri) from the dwarf locus for stars veiled in both the g and r bands is shown here 
against the g band veiling (r^). The models used are those shown in Figure 5. 
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r-i 

Fig. 7. — (r — i,r) color- magnitude diagram.. The Orion PMS locus is traced by the 
variable stars {squares: g band a > 0.05 magnitudes, diamonds: both g band a and z band 
a > 0.05 magnitudes.) in this color-magnitude diagram. The solid line indicates the main 
sequence at the m — M = 8.22 distance of Orion OBlb. The approximate survey bright 
limit based on saturation in the z band is marked by the dashed line and the distance to the 
background Gould Belt population is delimited by the two dotted lines. 
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Fig. 8. — {i—z, i) color-magnitude diagram.. The Orion PMS locus is traced in this color- 
magnitude diagram following Figure 7 above. The solid line indicates the main sequence at 
the m — M = 8.22 distance of Orion OB lb. The approximate survey bright limit based on 
saturation in the z band is marked by the dashed line and the distance to the background 
Gould Belt population is delimited by the two dotted lines. 
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RA 

Fig. 9. — Orion spatial distribution. The distribution in RA and DEC is shown here 
for the pre-main sequence candidates in the Orion sample. The variable stars are marked 
following Figure 7. A portion of the Orion OBlb assocation is visible for 82° < ci;2ooo < 85°. 
The NGC 2068/2071 protocluster is seen at 0:2000 = 86.5°. The cluster at 02000 ~ 81°, 
^2000 ~ 1° is part of the Orion OB la association. The lines trace the approximate boundaries 
of the portions of the L1630 cloud {upper left), the Orion OBla association [upper right), 
and the Orion OBlb association {lower half) included in this study. 




AQgri 

Fig. 10. — Shift in Qriz- In this figure we show the shift in Q^iz measured between the 
bright and faint states plotted against the J band magnitude. The dashed fines are the 3 a 
fimits based on the maximum photometric errors in the survey. Variable stars arc marked 
following Figure 7. From comparison with Figure 2 we see that a change in Q^z of 0.1 
corresponds to ~ 100 K in T^ff. 




AQgri 

Fig. 11. — Bright state AQg^. The bright state AQgri is shown here measured relative to 
the dwarf locus. These values correspond to g band veilings less than 1 (see Figure 6). The 
dashed line at AQgri = -0-7 indicates the color shift due to surface gravity effects alone for 
a pre-main sequence star having Tg// = 2700 K. Variable stars are marked following Figure 
7. The dashed lines centered around AQgri — are the 3 a limits based on the maximum 
photometric errors in the survey. Evidence for weak veiling is found in stars fainter the 
substellar boundary. 




Fig. 12. — Orion g band variability. The standard deviation of the g band variation 
is plotted against J for the Orion pre-main sequence sample. Variable stars are marked 
following Figure 7 above. The median ag values computed in one-magnitude bins are shown 
by the connected open triangles. The dotted line marks the HBL at 2 Myr for Orion and 
the dashed line traces the ag = 0.05 threshold for inclusion in this sample. The box in the 
lower right-band corner is defined by J > 14 and ag > 0.5. A lack of faint objects having 
high-amplitude variabilities is evident. 




CJg/az 

Fig. 13. — Ratio of g band to z band variability. The ag to az ratio is plotted against 
J for the Orion pre-main sequence sample. The median crg/cFz ratios computed in one- 
magnitude bins are shown by the connected open triangles. With the exception of 5 outhers 
in the lower right we see that the g band variation decreases relative to that in z for the 
fainter objects. The dotted hne marks the HBL at 2 Myr for Orion. The Ug — locus 
is shown by the vertical dashed line. The box in the lower right-band corner is defined by 
J > 14 and Ug/uz > 3. A relative lack of objects in this selection region may refiect a 
reduction in accretion shock temperature. 
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Fig. 14. — Orion (J — Kg, J) color- magnitude diagram. The 2MASS J magnitude is 
plotted against J — Kg color with the BCAH98 2 Myr and 20 Myr isochrones, converted 
from the CIT to the 2MASS system, shown for comparison. The dotted line marks the 
substellar boundary at this age and the diagonal line indicates the Ay = 5 reddening vector. 
The objects below the reddening vector are substellar candidates. Variable stars are marked 
following Figure 7. 



-36- 



Table 1. Orion OBI Sub- Associations 



Region Distance in pc (m-M) Age (Myr) 



Orion OBla 


330 (7.59) 


11.4 ± 1.9 


Orion OBlb 


440 ( 


8.22) 


1.7 ± 1.1 


Orion OBlc 


460 ( 


8.31) 


4.61-? 


Orion OBld^ 


450 ( 


8.27) 


< 1 



"The Orion Nebula Cluster is not included in this study. 



Table 2. SDSS Orion Imaging Runs^ 



Run 


Date (MJD) Interval (Days) Elapsed (Days) 


0211 


51115 




0259 


51134 


19 19 


0273 


51136 


2 21 


0297 


51139 


3 24 


0307 


51140 


1 25 


2955 


52312 


1172 1197 


2960 


52313 


1 1198 


2968 


52314 


1 1199 


4158 


52912 


598 1787 


'^This table lists the SDSS imaging of the Orion equatorial region giving the run number and the MJD of the observation. 


The Interval and Elapsed columns indicate the temporal sampling by giving the number of days between successive observations 


and the total number of days since 


the first imaging, respectively. 




Table 3, 


. Limiting Magnitudes and Photometric Errors^ 


Band 


Survey Faint Limit 


Maximum Error 


9 


SDSS 21.5 


0.05 


r 


SDSS 20.0 


0.03 


i 


SDSS 18.0 


0.02 


z 


SDSS 17.0 


0.02 


J 


2MASS 15.5 


0.05 


H 


2MASS 14.5 


0.05 


Ks 


2MASS 14.5 


0.10 



"These are the faint magnitude limits and maximum photometric errors observed when the selection criteria described in the 
text are applied to the SDSS and 2MASS data acquired in the Orion equatorial region. 



Table 4. Scaling Between Actual and Observed a ^ 



N Epoch Factor (99%) 



2 


1.748 


3 


1.344 


4 


1.226 


5 


1.155 


10 


1.072 


15 


1.040 


20 


1.033 


25 


1.026 


30 


1.021 



''This tabic was created, as described in the text, by a series of Monte Carlo simulations to determine the maximum standard 
deviation (a) in a random sequence given an intrinsic a of the parent population and a limited number of observations (N 
Epoch). The scaling factors reflect the maximum a seen in 99% of the simulated data sets at each value of N Epoch. In the 
limit of large N Epoch the scaling factors tend towards unity. 



